Existence of Spherical Initial Data with Unit 
Mass, Zero Energy, and Virial less than - 1/2 
for the Relativistic Vlasov-Poisson Equation 
with Attractive Coupling 



In a recent paper, Kiessling and Tahvildar-Zadeh proved that any 
classical solution of the relativistic Vlasov-Poisson equation with at- 
tractive coupling launched by spherically symmetric initial data with 
zero total energy and virial less than or equal to -1/2 will blow up in 
finite time. They left open whether such data exist. Subsequently, the 
question was raised whether any such data exist at all. In fact, the 
simplest conceivable ansatz, a nearly uniform ball of material centered 
at the origin with momenta directed inward, must have virial strictly 
larger than -1/2! In this brief note, we settle this issue by constructing 
two classes of such initial data. 



1 Introduction 

The relativistic Vlasov-Poisson (rVP) system is given by 
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rVP + models a system with repulsive interaction while rVP~ models a system 
with attractive interaction. One of the earliest papers to appear on the 
subject is [GS85] wherein Glassey and Schaeffer show that global classical 
solutions to rVP ± will exist for initial data that are spherically symmetric, 
compactly supported in momentum space, vanish on characteristics with 
vanishing angular momentum, and have £°°-norm below a critical constant 
C± with C+ = oo and < oo. Recently, Kiessling and Tahvildar-Zadeh 
[KTZ08] have extended the theorem of Glassey and Schaeffer for rVP~ by 
proving global existence of classical solutions for the same type of initial data 
except that they only need to have ii^-norm below a critical constant with 
Cp < oo and identically zero for < 3/2. The sharp value of this constant 
was subsequently determined for all values of (3 in [You 10] . 

Glassey and Schaeffer also investigated what may happen when rVP~ is 
launched by initial data with ||/||oo > C~ . They proved that negative energy 
data lead to "blow-up" in finite time. In the same paper extending the global 
existence result, Kiessling and Tahvildar-Zadeh proved that any spherically 
symmetric classical solution of rVP - launched by initial data /q G ^3 fl C 1 
with zero total energy and total virial less than or equal to -1/2 will blow up 
in finite timeS. Specifically, with the total energy of fo given by 

£(/o) = JVl + WfiP, -\JJ f ^\™4\ q) d ^'i (!) 
and the virial of fo given by 

V(/o) = f(q-p)fo(p,qW, (2) 



where J ■ ■ -dfi = J J ■ ■ -d 3 pd 3 q, Theorem 6.1 of [KTZ08] states: 

Theorem 1. Let t — > f t G ^3 fl be a spherically symmetric classical 
solution of rVP~ over some interval [0, T) launched by initial data fo with 

3 /15^ /3 



||/o||3/ 2 > g ( 16 



IfS(fo) = and V(/ ) < -1/2, then T < 00 



^By CP3 n C 1 we mean the set of probability measures on K 6 whose first three moments 
are finite, that are absolutely continuous with respect to Lebesgue measure, and whose 
Radon-Nikodym derivative is £} . 
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The proof is by contradiction. The zero-energy condition ensures that the 
virial is strictly decreasing. Assuming that an initial datum satisfying the 
requirements above launches a global solution, one can show (using a number 
of technical estimates) that a strictly positive quantity is bounded above by a 
sum of positive terms and a time integral of the virial. The authors then show 
that after a sufficiently long time this strictly positive quantity is bounded 
above by something negative. Hence, such an initial datum cannot lead to 
global existence in time. 

Now, the authors of |KTZ08] did not give any specific examples of fa- 
vorable initial data, and recently the question was raised whether such data 
actually exist |CallO] . Indeed, nearly uniform balls of material centered at 
the origin (the simplest ansatz one could conceive) cannot be made to work. 
Surprisingly, the existence of these initial data is an interesting, non-trivial 
question! The purpose of this brief note is to settle this issue. 

We first define a generic separation of variables ansatz and derive all basic 
formulae needed to carry out our subsequent investigations. We then show 
that the simplest ansatz (a nearly uniform ball of material centered at the 
origin) cannot have the desired properties. Next, we employ a disjoint core- 
halo structure (a central ball of material with a separate thin ring of material 
further out) to construct a class of initial data satisfying the requirements in 
[KTZ08J. In addition, we show that zero-energy initial data can have arbi- 
trarily negative virial. The success of the core-halo scheme raises the question 
of whether such a disjoint structure is necessary for our requirements. We 
answer this in the negative by exhibiting a second class of initial data with 
zero total energy and total virial less than or equal to -1/2 that have the nice 
properties of being monotonically decreasing and "singly supported." 

We should mention that the details of the blow-up scenario are not ad- 
dressed in either [GS85J or [KTZ08J. However, in |LMR08j Lemou, Mehats, 
and Raphael proved that "blow-up" for systems launched by initial data with 
negative total energy approaches a self-similar collapse. 

Lastly, we also mention other recent results on rVP~. Hadzic and Rein 
[HR07] showed the non-linear stability of a wide family of steady-state solu- 
tions of rVP~ against a class of perturbations utilizing energy- Casimir func- 
tionals. Shortly thereafter, Lemou, Mehats, and Raphael II.M HOSl ILMR 09J 
investigated non-linear stability versus the formation of singularities in rVP~ 
through concentration compactness techniques. While these results are inter- 
esting in their own right, they will play no role in the following considerations. 
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2 The Separation of Variables Ansatz 

We first establish the various formulae we shall need for our particular class 
of functions. Consider a generic separation-of-variables ansatz: 

f(p,q) = C V (\q\M\p\)J?(cos(9 p , q )), (3) 

where 6 p%q is the angle between q and p (both considered as vectors in M 3 ) and 
77, and ££ are non-negative. The coefficient C will be chosen to normalize 
the total mass to 1. 

The total mass is given by 

/•OO /*00 1*1* 

f( P ,qW=c v(\q\)\q\ 2 d\q\ m P \)\ P \ 2 d\p\ // j2?(cos(0 M ))dnpdn, 

Jo Jo JJ&xn 2 

/'CO poo pi 

= 8n 2 C / v(\q\)\q\ 2 d\q\ / ®(\p\)\p\ 2 d\p\ / 5?{x)dx. 
Jo Jo J-i 

Using the notation 

POO 

\\g\\i = / \9(r)\dr, 



and recalling that all functions under consideration are non-negative we see 
that taking 



C- 1 = 87r 2 \\ V \q\ 2 \\ 1 \\<5>\p\ 2 \\ 1 j 



J>?(x)dx 



will ensure that / has total mass 1. Similarly, the £ 3 / 2 -norm is given by 

2/3 



3/2 



JT v(\q\) 3/2 \q\ 2 d\q\ Jo" H\p\) 3/2 \p\ 2 d\p\ /i ^{xf l2 dx 
27r 2 /3|| r? |g|2|| 1 ||$| J9 |2|| 1 j^^( a .)^ 



Recalling the definition of the virial functional ([2]) from the introduction, 
for our separation of variables gives 



V(f) 



poo poo r p 

C V (\q\)\q\ 3 d\q\ m P \)\p\ 3 d\p\ // cos(6 p , q )J?(cos(9 p>q ))dn p dn q 

J0 J0 JJs 2 xS 2 



87rX||r7|g| 3 ||i||$|pH|i / x££{x)dx 

\\v 



-1 

'x 



q\ 3 h \\^\p\ 3 \\i!\ x ^( x ) d ^ , 
M?l 2 lli||$|p| 2 lli J 1 _ 1 ^(x)dx ' 1 j 
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Of course, in order to have a negative virial we must require 



J xJf(x)dx < 0. 



For the energy (as given by ([I])), we compute the "kinetic" and poten- 
tial energy contributions separately (the portion of the energy we label as 
"kinetic" also contains the rest mass energy). For the kinetic energy, we have 

KE(f) = J VT+Wf(p,qW 

roc poo pp 

= C H\p\)y/T+W\p\ 2 d\p\ h(\q\)\q\ 2 d\q\ h?(cos(6 m ))dn p dn q 

JO JO JJs 2 xS 2 

_ \\*VT+W\pFh m (5) 



||$b| 2 ||i 

We begin computing the potential energy by first computing the spatial 
distribution p associated to /: 

p(q) ee J f(p,q) d 3 p 

oo 



Cr,(\q\)J H\ P \)\p\ 2 d\p\ I J*f(cos(0 M ))dfi f 

v(\q\) 



o 



47r||?7|g| 2 ||i 

In terms of this function, the potential energy contribution is now given by 



\q - q' 



We calculate this integral directly, and after a little work find 
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PEU) = -JTum ^M)M(/ V(W\)W\ d\q \ ) d\q\. (6) 
\\V\Q\ Hi Jo 
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3 Nearly Uniform Balls Fail (Barely) 



At this point, we can show that the simplest ansatz cannot work. Relaxing 
the differentiability requirement for a moment, we make the specific choices 

v(\q\) = X[o,R](|g|), 

H\p\) = X[o,p](bl), 

= X[-l,a]( X ) 

where R > 0, P > 0, — 1 < a < 1, and xi is the characteristic function of the 
interval I. We have that 



KE{f) 
PE(f) 



The zero-energy requirement allows us to solve for R in terms of P. Cal- 
culating the virial and plugging in the formula for R = R(P) gives V(/) in 
terms of the parameters P and a. Simple asymptotics shows that for this 
ansatz 

v(/) > -I 

for any choice of parameters. Smoothing out the boundaries of these step 
functions (taking care to keep the total energy zero) shows that this ansatz is 
untenable. Of course, this type of initial condition may well lead to collapse 
after finite time, but this cannot be verified by our virial condition. 

The failure of nearly uniform balls to satisfy our criteria is a little sur- 
prising at first, but upon reflection we can understand the difficulty Recall 
that the "kinetic" energy term contains the rest mass contribution and so is 
at least one (in our rationalized units). Heuristically, this forces the radius 
of the ball of mass to be small in order to have the potential energy negative 
enough to balance the "kinetic" energy contribution. This has severe conse- 
quences for the virial - which is made smaller in magnitude by concentrating 
all the matter near the origin. 
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4 Success: A Disjoint Core-Halo Ansatz 



The considerations at the end of the previous section lead us to try a disjoint 
core-halo scheme: 

V(\q\) = X[o,iii](|?|) + aX[ife,Jfe](l?l), (7) 

H\p\) = X[o,p](bl), (8) 

2{x) = X [-i,a]( x ) (9) 

where < Pi < P 2 < P3, < a, < P, and — 1 < a < 1. Again, we have 
relaxed the differentiability requirement to make the computations tractable. 
The idea behind this ansatz is that putting a tiny halo of material far from the 
center should increase the magnitude of the virial while keeping the potential 
energy negative enough to maintain zero total energy for the system. 
We first consider the specific choices 

Pi = -, P 2 = 1, P 3 = 2, and P = 1. 
5 

The zero energy condition forces (thanks to Maple) 

35 ln(l + y/2) + 30 - 105^ + 2^6480^ - 1655 - 2160 ln(l + y/2) 

125 (735 - 188 - 245 ln(l + y/2)) 

(a small but positive number). Again by Maple, choosing — 1 < a < —4/5 
gives V(f) < —1/2. Since our ansatz is compactly supported, we have the 
requisite number of moments. To complete the argument, we note that we 
can smooth out the various step functions in such a way that the resulting 
integrals are as close to the values obtained above as we like. Since a was 
chosen to force the zero-energy condition and the necessary value was strictly 
positive, it follows that by smoothing the test functions appropriately we can 
still choose an a > to keep the total energy zero. Using Maple to check 
the £ 3,/2 -norm for the range of a listed above shows that our data are well 

1 /3 

above the critical norm | (ff) • Thus, we have constructed an initial datum 
satisfying all the requirements of Theorem 6.1 in [KTZ08j. 

Considering the lower bound on the virial for nearly uniform balls, the 
reader may wonder whether a similar phenomenon holds for the disjoint core- 
halo structure. We give an argument showing that the virial can be made 
arbitrarily negative in this case. We make the following choices in ([7]): 

Pi = P" 2 , P 2 = P, and P 3 = P 2 . 



7 



Now, the zero energy condition allows us to solve for a = a(P). Using Maple 
for convenience, one sees that as a function of P, a is positive for sufficiently 
large P and is proportional to p~ 23 / 2 as P — > oo. Plugging these choices into 
the formula for the virial and looking at the asymptotics for large P shows 
that the virial is proportional to — (1 — a)P 3 . Smoothing these functions out 
as above shows that the virial is unbounded below. As before, the £ 3//2 -norm 
for these functions remains well above the critical norm for a wide range of 
values of the parameter a. 



5 A Monotonic Core-Halo Ansatz 

Having established the existence of initial data satisfying the theorem of 
Kiessling and Tahvildar-Zadeh, we are technically done. However, we can 
ask whether the separation of the halo from the core is an absolute necessity. 
The answer is no. Even though nearly uniform balls (after regularization) 
were shown not to work, we can find zero-energy initial data with virial 
less than —1/2 that are monotonically decreasing and supported in a ball 
centered at the origin. 

As in the previous section, we describe initial data that are not smooth 
with the understanding that we can appropriately regularize our functions 
without significantly altering any of the required properties. We now consider 
the ansatz 



= X[o,Rx](\<i\) + f^yj X[flx,jfe](M) + Qr] X[r 2 ,r 3 ](\<i\)i ( 10 ) 
*(bl) = x[ ,p](H), (ii) 

J2f(x) = X[-i, a] ( x )- (12) 

Note that we have kept the same choices for the momentum and angular 
portions of / but substituted a spatial portion that is monotonically decreas- 
ing inside the ball of radius R3 centered at the origin and zero outside this 
ball. A representative diagram for 7] is given in Figured] below. Note that it 
retains the essential ingredient from the previous section. Namely, there is 
a dense ball of material centered about the origin (to sufficiently lower the 
potential energy) surrounded by a thin atmosphere of material further out 
(to sufficiently increase the magnitude of the virial). 

Even with this relatively simple ansatz, the calculations are unwieldy. 
Thanks again to Maple, we merely report a choice of parameters that does 
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Figure 1: Representative monotonically decreasing r\ 



the job. Namely, choosing 

i?i = , Ro = — i = — i ti = 3 

100 11 10 

sets a corresponding P (about 19.69) to force the zero-energy condition. 
Finally, any choice of a less than roughly —9/10 gives us a virial less than 
-1/2. A check of the £ 3 / 2 -norm shows that it is well above the critical value. 
We remark that our choice of parameters was totally arbitrary, and so we 
expect that an entire range of such parameters will do the job. However, 
here we are content merely with the existence of one such solution! 
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